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■ The ground state single particle Green's function describing hole propagation is calcu- 

lated exactly for the 1/r 2 quantum many body system at integer coupling. The result is in 
agreement with a recent conjecture of Haldane. 
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1 Introduction 

The 1/r 2 quantum many body problem, and the closely related Haldane-Shastry 1/r 2 Heisenberg 
chain, are the subject of intense theoretical study at present, due in part to their relationship to 
an ideal gas obeying fractional statistics (see [1] for a review). One direction of study has been 
the exact calculation of some ground state correlation functions [1-5] . Let us briefly summarize 
O ■ the main results of these works. 

. The static ground state correlations (one-body density matrix, two-point distribution func- 

tion) for the 1/r 2 Hamiltonian in periodic boundary conditions 

gg + ^-'Hi) J^ STK^P: -' 6 ** (L1) 



(here we have set h 2 /2m := 1), for which the exact ground state wave function is 

\0 >:=ip (xi,. . . ,x N ) := J[ (simr(x k -Xj)/Ly, (1.2) 

l<j<k<N 

have been calculated exactly in terms of a class of generalized hypergeometric functions in 
several variables based on the Jack symmetric polynomials [2]. These functions are denned in 
terms of power series, and in the cases of interest an explicit expression for each coefficient is 
known. However, the coefficient of the iVth term involves a sum over all partitions of N into q 
parts and is thus impractical to compute except for small N and q. To overcome this difficulty, 
recurrence formulas have been presented [3] which allows for the rapid numerical evaluation 
of the ground state correlations in the finite system. Furthermore, integral representations of 
the generalized hypergeometric functions are known [4,5], from which the density matrix and 
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two-point distribution function have been written in terms of q and 2q dimensional integrals 
respectively. 

The thermodynamic limit can be taken in the integral formulas, giving similar integral for- 
mulas for the density matrix and two-point distribution function in the thermodynamic limit 
[5] . This integral formula for the density matrix has recently been generalized by Haldane [1] to 
an integral formula for the retarded single-particle Green's function 

- iG({x, t), (0, 0)) := TL < | tp + (x, i)^(0, 0) | > (1.3) 

where TL denotes the thermodynamic limit N,L — > oo, N/L = p. Haldane conjectures that for 
general q € Z + 

-iG((M),(0,0)) =\B q J d Vl ...dv q f[(l-v])- 1+1 '* 



-1.1] 



j=i 



x H \v k - Vj \ 2 /« J| e »w e -^p) 2 (i-^ 2 )^ ( L4 ) 
i<i<fc<<? i=i 



where 



As t — > + the results of [5] for the density matrix are reclaimed. For general t, this result was 
proved in [6] for q = 2 using a mapping to a dynamical matrix model; for q = 1 - free fermions 
- it is easy to verify from the differential equation satisfied by — iG^(x,t), (0,0)^ (see e.g. [7]). 

2 Derivation of Haldane's Conjecture 
2.1 Action of the Hamiltonian on a hole state 

Our starting point is to write the retarded single-particle Green's function for the finite system 
in a more explicit form (see [6], eq. (12): 



N 

iG N+1 ((x,t); (0,0)) = (N + l) < | H(smir( Xl -x)/L) q e -^-E )t/n 

i=i 



N 



x Y[(smirxi/L) q \ > (2.1) 
l=i 

Next we introduce the auxiliary variables y±, . . . ,y q and take up the task of calculating 

N q N 

(N + 1) < | J] Y[ smTr( Xl - yj )/Le~ i( - H ~ E ^ h J^sinvr^/L) 9 | > . (2.2) 
l=i j=i i=i 

We do this by first expanding the exponential 

e -i(H-E )t/H = g H*/E) J {H _ Eq) j (2 3) 

3=0 J ' 
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and considering the action of the operator H — Eq on the hole state 



\(j) >:= Yl sin7r(xi - yj)/L | > . 



(2.4) 



We have the following result. 
Lemma 1 

With H given by (1.1), Eq the corresponding ground state energy and | > given by (2.4) we 
have 

(H - Eq) I <j> > = T {y} | (j) > (2.5a) 

where 

' « d 2 IT 2 « 



T {y} : = 1 

( d 



+ 



7T 



(qN + N(N - 1)) 



(2.56) 



Proof 

Direct differentiation gives 

(H-E )\(f>> 

2 



N q 

qN + Y^ cot7T (xi - y h )/L cot tt{xi - y h )/L 



+ 2 lYl cot7T (xi - x v )/L ^cot7r(x/ - yj)/L 

3=1 



=1 f'=i 



> 



(2.6) 



We rewrite the second summations by grouping together the summand with the summand with 
I and I' interchanged. Simple manipulation gives 

cot 7r(xi — xii)/L [cot tt(xi — yj)/L — cot tt(xi> — yj)/L] 

= - (1 + cot7r(x/ - yj)/L cot ir(xi> - yj)/L) . (2.7) 
In the first summations we use (2.7) to rewrite the summand: 

cotvr(^ - y h )/L cot ir(xi -y h )/L 

= -1 - cotvr^ - y j2 )/L[cotTr(y jl - x t )/L - cot7r(y i2 - x t )/L\. (2.8) 
Substituting (2.8) and (2.7) in the first and second summations respectively of (2.6) we obtain 

{H - Eq) I <f> >= - ( - ) [ - q 2 N - <***{Vji - y j2 )/L 



N 



Y ( cot7T (yji - X l)/ L - cot7T (y j2 - Xi)/L) - qN(N - 1) 



1=1 



q N 



-qJ2Yl cot7 r(xi - yj)/L cot tt(xi> - yj)/L] 



(2.9) 
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But from (2.4) and (1.2) 



- 10 >= - (Jj cot < x i - yi)i L 1 > ( 2 - 10 ) 

I > • (2.11) 

□ 

Remark: The operator T^ y y has occured previously in the study of the 1/r 2 quantum many 
body system [8]. Thus if ip denotes an eigenfunction of H with q replaced by 1/q (recall (1)), 
and ijj = $ | >, then $ is an eigenfunction of T y . Furthermore, it is known [4] that the 
eigenfunctions of T„ are the Jack polynomials c[ q) (e 2niyi / L , e 27tiy "' L ) where k is a partition 
which labels the eigenfunction. This latter fact will play a crucial role in our subsequent analysis. 



and 



d 2 , . (^ 2 N 



.. , i r >= ( y XI cot7r ( x i - Vj)/ L cotvr(x// 
KLJ i', i=i 

Substituting (2.10) in (2.9) we obtain (2.5), as required. 



Using Lemma 1 and (2.3), we see that 

-iG N+1 ( (x,t),(0,0) ) 



(#+i)E 



( i T ^(T M ) J <<AlIl sin W^k> 



yi=-=y q -- 



(2.12) 



We have previously evaluated the inner product in (2.11), which at the specified point is precisely 
the static one-body density matrix. Thus [4, eq.(3.4)[J 



N 



(N+ 1) < <f> | Y[(smTTXi/L) q | >= 



N + 1 



C,,jv II e'^y^Fi^i-N, 1--N + 1- 1/q; e 2niyi / L , . . . , e 2niy "/ L ), (2.13) 



where 



with 



a 



q,N 



A q (l/q,-l,2/q) 
AJl/q,-N- 1,2/q) 



(2.14) 



r(l + A/2)r(Ai + A 2 + A(n + j - 2)/2) 



nl 2 ' j /i r(i + Aj/2)r(A 1 + A(j - i)/2)r(A 2 + a(j - 1)/2) • 

The generalized hypergeometric function of several variables 

(a, b; c;xi, . . . , x m ) is the unique symmetric power series solution of the partial differential 
equations 

d 2 F r 2, _ . 2. ,9F 

-(m- 1) - (a + 6+ 1 - -( 
A A 



+ 



E 



i 



fc=1 £fc a-j 



d 

xj(i - x j)-q^: - x k{t - X k ) 



_d_ 

dxi 



F-abF = 



(2.15) 



1 this equation is erroneously missing the factor C q ,N 
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with the initial condition 2-^1^ = 1 when x\ = . . . = x m = [9,10] (when m = 1 - single variable 
case - the differential equation satisfied by the Gauss hypergeometric equation results). 

We need to be able to compute the action of the operator (T^y on (2.12). For this purpose 
we first write 2-^1^ in integral form. 

2.2 An integral formula for 2^1^ 

The following integral representation is due to Zan [10]: 

2 F 1 W x \a, Ai + A(n - l)/2; Ai + A 2 + A(n - 1); *i, . . . , z n ) 

= A n (X 1 ,X 2 ,X) / i^" (2/A) (a;^i,---,^n;si,---,Sn) 

Jin 

x D Xu x 2 ,\(si, s n )dsi ...ds n (2.16) 



where 



and 



Dx lM ,x{8 U ...,8 n ):=n*i(X-*i) X *~ 1 II \*k-s 3 \ X (2-17) 

1=1 l<j<k<q 



i^o (2/A) (a; 



Z\ j ■ ■ ■ > Z n ! ^1 j 



00 ^(2/A)/ \^2/A/ \ 

E \ " r„l(2/A) °^ \Zl>---: Z n )^K [Si,..., S n ) fn tq\ 



fc=0 ■ |«|=fc ^ (,1, 



In (2.18) the second sum is over all partitions («i, . . . , re n ) of A; into n parts, cl 2 ^ denotes a 
suitably normalized Jack polynomial [9], and 



a 



(a) - 



1 



ni --^- 1 )) ( 2 - i9a ) 



where 

(a) fc :=a(a + l)...(o + fc-l). (2.196) 

In general no formulas expressing ±J-q ( 2//- ^ in terms of simpler functions are known. However, 
an exception is the equal variable case z\ = . . . = z n = z when [4] 

n 

!^ 2/A) (a;z,...,z: Sl ,...,s n ) = l[(l- ( 2 - 2 0) 

The interval of integration for each variable si,...,s n given in [10] is / = [0, 1]. 
The l.h.s. of (2.16) corresponds to the 2 -Fi ^ function in (2.12) if we take 

2/X = q, a = -N, n = q, X 1 = l/q, A 2 = -N - 1, Xi = e 2niyi/L (2.21) 

in the former. However, we see from the definition (2.17) that in this case the integrand in (2.16) 
is not integrable at Sj due to the factor (l—Vj)~( N+2 \ for each j = 1, . . . ,q. In this circumstance, 
in the one variable case (Gauss hypergeometric function) we know that the contour of integration 
/ can be shifted to 

J = (-oo,0] (2.22) 

and the integration formula (2.16) is both well defined and satisfies (2.15). We expect this 
property to carry over to the m- variable case, and indeed it does with the change of normalization 

A n (X 1 , A 2 , A) -> (-l) Al M n (A!, -A(n - 1) - A x - A 2 + 1, A) (2.23) 

as is shown in Appendix A. Thus we can use the integral representation (2.16) with parameters 
given by (2.21), interval of integration (2.22) and change of normalization (2.23). 
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2.3 An equivalent operator 

Using the integral representation in (2.12) we see from (2.16) and (2.11) that we must calculate 



l=i 



{ T {v}) 

To accomplish this task we recall [4] that T^ y y has 

i 



yx=...y q =x 



(2.24) 



1=1 



as an eigenfunction, with eigenvalue t K say, for each partition k. Hence, using (2.18), we see 
that (2.24) is equal to 



oo -, r <(q)(„2mx/L Jlmx/L^W) / \ 

e -7Ti P qX ±_ [g] (2/X) (t K y ^ '•••' e ' )CK {S 1 ,...,Sg) 



3=0 h \k\=j 



d q \i,...,i) 



(2.25) 



If we know introduce an operator T| s j acting on the coordinates s±, . . . ,s q such that 

T[ s} C^( Sl , ...,s q ) = t K C<s\ 8l , ...,s q ) (2.26) 

we see from (2.18) that 

(T[ s] y e-^^i-N; e 2mx > L , e 2mx / L ; s u . . . , s q ) (2.27) 



is identical to (2.25), so we can replace (2.24) by the equivalent operation (2.27). The con- 
struction of Tj s | is simple, due to the symmetry between the eigenfunctions cj^\z±, . . . , z q ) and 



Ci 9) (si, . . . , s q ) exhibited in (2.18). In Appendix B we show that 

(t) 



+ (q(4n 2 p/L) - ( 2 vr/L) 2 ) £ s~ + q(l 

j=i i 



(2.28) 



The calculation of (2.27) is straightforward since the special case of the function therein 

can be evaluated according to (2.20). 

Substituting the evaluation (2.20) in (2.27) and replacing (2.24) by the resulting formula we 
thus have from (2.11), (2.12) and (2.16) (with the further specifications made after (2.23)) the 
formula 

-iG N+1 ((x,t),(0,0)) 

-it/ny 



-C N , q A g (l/q, N + l/q, 2/q)e-* i "> x £ ■ ., 

j=0 J - 



L 



x / dsi . . . ds q 

'(-OO,0]9 



{Tls } ) 3 I[(l-e 2mX/L si) N 



i=i 



X D l/q-N-l,2/q( s li ■ ■ ■ , s q) 



(2.29) 
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2.4 The thermodynamic limit 

From the explicit form (2.28), the action of the operator (t^ s ^J' > in (2.29) can readily be com- 
puted using a computer algebra package for given N and j. However, many terms result and it 
doesn't seem possible to sum the series in j. However, an analytic calculation is possible in the 
large N, L limit. From (2.28) we see that to leading order in N with respect to its action on the 
function 



TJ(1 _ e 2 ™/ L Sl ) N (2.30) 
l=i 

we can replace (^{ s }) by 

The operator (2.31) gives terms which are 0(1) and 0(1/N) whereas the remaining part of 

j 



T| s |J gives terms 0(1/ N) only, which can therefore be ignored. The 0(1) action of (2.31) on 
(2.30) is readily computed, and we conclude 



{T' {s} yw i -^ ix,Ls i) N 
i=i 



^P) 2 E TTTTH + -rrr) IK 1 " e 2 ^ L Sl ) N + 0(1/N) (2.32) 

3=1 ^ S 3> 1 S 3 j l=i 

Substituting (2.32) in (2.29) we recognise that the sum over j is simply the power series for the 
exponential. Hence 

-iG N+1 ((x,t),(0,0)) 

~ -C N , q A q (l/q,N +l/q,2/q)e-^ x 



f q 
X / ds\ . . . dsq TT exp 

y(-oo,o]« " 



£ _ e 2 ™l L Sl ) N 



(i- Sj y i-sj 

x D 1/q _ N _ lt 2/ q (si,...,s q ) (2.33) 

We are now close to deriving the result (1.4). 

To finish off the calculation, we first note from (2.13) and (2.14) that 

lim C N , q A q (l/q, N + 1/q, 2/q) = B q (2.34) 

N— »oo 

where B q is given by (1.4b). Next we change variables 

= -tj (2.35) 

in the integral (2.33) so that it reads 

r q r i 

- / dh . . . dtg TT exp it/h(t 2 - tj) (1 - tj(l - e 2 ™/ L )) N 

J [0A] q i=i 1 1 
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x D l/q,l/q,2/q(h, ■ ■ ■ ,t q ) 



(2.36) 



But 



N,L — >oo 
N/L=p 



lim (1 — tj(l — e 




(2.37) 



Substituting (2.37) in (2.36), changing variables 



h = (vj + l)/2 



(2.38) 



and substituting the resulting expression in (2.33), we obtain the formula (1.4) for the single- 
particle Green's function (1.3). 

Note added: Since completing this work, the generalization of the formula (1.3) to all rational 
values of the coupling q, together with a derivation using Jack polynomials, has been announced 
by Ha [12]. 
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Appendix A 



Here we derive the integral representation (2.16) with interval of integration (2.22) and change 
of normalization (2.23). We start with the transformation formula [4, eq. (3.6)^[] 

2 F 1 W x \a,b;c;t 1 ,...,t m ) 

= C 2 F 1 W(a, b;a + b+l + X(n- l)/2 - c; 1 - t h . . . , 1 - t m ) (Ala) 

where 

C:=CU(.,^A) := /f - X{n -T' C : b - X{ "- (A«) 
A n (o — A(n — l)/2,c — b — a — A(n — lJ/2, AJ 

which is valid for a € Z<o- Use of (Al) in (2.16) gives 

2 F 1 W X '> (a, Ai + A(n - l)/2; a + 1 - A 2 ; z 1} . . . , z n ) 

A„(Ai,A2,A) f (2/x) i -i -i \ 

= 7T7 — x , w lWo x , . , w TVTT / ^ >[a;l-zi,...,l-z n ;si,...,s n ) 

C n (a, Ai + A(n - lJ/2, Ai + A 2 + A(n - 1); A) ./[o,i]» 

x £>Ai,a 2 ,a(si, • • • , s„)dsi . . . ds n (A2) 

where again it is assumed a £ Z<q. 
Next we change variables 

Si = -- — — (i = l,...,n). (A3) 
1 -Ui 

To do this we note 

-Dai,a 2 ,a(si, ■ ■ ■ , s n )dsi ...ds n 
= (-l) Xrn D Xl -[A 1+ A 2 +A(n-l)-l],A(^l, • • , u n )dux . . . du n (A4) 
Furthermore we can make use of the transformation formula given by the following result 

Lemma 
We have 

i^b 1 1 '{a; 1 - z 1} . . . ,1 - z n ;- 



1-ui l-u n 



11(1 — % ) a i^0 (2/A) («; Zl, .... ... ,«n) (A5) 



Proof 

We know that [10] 

n 



2 Fi( 2 / A )(a,6; C ;y 1 ,...,yJ = n(l-%')~ a 2^i (2/A) Kc-6;c;--^,...,--^) (A6) 

/ = i l — i/i 



ables s,- i— > 1 — s,- on the r.h.s. shows 



Expressing both sides of (A6) in terms of the integral representation (2.16) and changing vari- 
es 

i^" (2/A) («; yi, ■ ■ ■ ,y n ', si, ■ ■ ■ , s n ) 

[0,1]' 



2 this formula is erroneously missing C 
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- II (1 - W 2/A) (a; -j^-, -j^-; l- Sl ,...,l-s n ) 

x -Dai,a 2 ,a(si, • • • , «n)dsi ...ds n = (A7) 

Since the terms in the square brackets are independent of Ai and A2 while (A7) is identically 
zero for all Re(Ai), Re(A2) > 0, we can conclude the combination of terms in the square brackets 
vanishes identically. The identity (A6) then follows by noting from (2.18) that in general 

i^" (2/A) (a;yi, ■■■,y n ;si,...,s n ) 

is unchanged by interchanging all the variables y±, . . . , y n with the variables si,...,s n . 
Substituting (A4) and (A5) in (A2), noting that in general 

n 

XT (1 -Uj) a Dx 1 ,\ 2 ,\{u 1 , ...,%) = Dx u \ 2 +a,\(ui, • • -,Un), (A8) 

J'=l 

and replacing A2 by a + 1 — (Ai + A2 + A(n — 1)) shows that (A2) reduces to (2.16) with the 
substitutions (2.22) and (2.23), as required. 

Appendix B 

Here we construct the operator T| s j defined in subsection 2.3 by the eigenvalue equation 

T' {s} Cj?\ Sl , ...,s q ) = t K Cj?\ Sl , s q ), (Bl) 
where the eigenvalue t K is first to be calculated from the eigenvalue equation 

T {y} l[ e -™pyi C M(e 2niyi/L , e 2niy * /L ) 
1=1 

Q 

= t K J] e- wipyi C^{e 27Tiyi / L , . . . , e 2niy o /L ) (B2) 
1=1 

with Tjy} given by (2.5b). 

To calculate t K we recall that the Jack symmetric polynomial 
ci?\e 2myi / L , . . . , e 2myc i 1 ' L ) , which is labelled by a partition k = (k±, . . . , n m ) of non-negative 
integers such that 

m 

k\ > K2 > ■ ■ . > K m and ^2 K j = (B3) 

3=1 

is the unique (up to normalization) solution of the eigenvalue equation [4;eq. (2.11) with a = 
m = q] 

T {y} Ci g \e 2niyi / L ,...,e 2niy "/ L ) 
= (^qeUq) + {jf[qN + N(N-l)]j C^(e 2my ^ L , . . . , e 2 ^/ L ) (B4) 

where 

e'M = -( 2 -^]\eM + k/2) (B5a) 



q \LJ 

eM = q E Kj(Kj - l)/2 - ]T(j - l) Kj + (q- l)k (B56) 

3=1 q=l 
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Using (B4) and the fact that is homogeneous of order k, which implies 

q n 

V C ( K q) (c 27Tim/L , e 27r ^ /L ) 
3=1 d0j 

= ^l kC ^( e 2 ^ L , e 2m y^ L ), (B6) 

Li 

the l.h.s. of (B2) is readily computed, and we find 

t K = 2(27r/L) 2 e K (q) + k[q{^ 2 p/L) - (2tt/L) 2 } + q(q - 1)(tt/L) 2 N (B7) 
To construct T^ g y with this eigenvalue acording to (Bl), we recall [9] that 

A w C(«>(*i, ...,s q )= e K (q)C^( Sl , . . . , s q ), (B8) 

where 

<7 <-i2 <? -? 2 9 

Furthermore, we can rewrite (B6) as 

£ s~C^( Sl , ...,s q )= kCj?\ 8l , ...,s q ) (BIO) 

j=l 3 

Hence 

( - 2 (2vr/L) 2 A {s} + [g(4^ 2 p/L) - (2vr/L) 2 ] ]T ^ 

j=i J 

+g( (7 -l)(7r/L) 2 iv)cW( Sl ,..., Sg ) 

= (Bll) 
from which we read off that the explicit form of T^ s y is given by (2.28). 
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